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A thermodynamic model is presented to describe possible mechanisms of ice formation 
on unheated surfaces inside a turbofan engine compression system from fully glaciated ice 
crystal clouds often formed at high altitude near deep convective weather systems. It is 
shown from the analysis that generally there could be two distinct types of ice formation: (1) 
when the ‘"'surface freezing fraction’^ is in the range of 0 to 1, dominated by the freezing of 
water melt from fully or partially melted ice crystals, the ice structure is formed from accre- 
tion with strong adhesion to the surface, and (2) when the “surface melting fraction” is the 
range of 0 to 1, dominated by the further melting of ice crystals, the ice structure is formed 
from accumulation of un-melted ice crystals with relatively weak bonding to the surface. 
The model captures important qualitative trends of the fundamental ice-crystal icing phe- 
nomenon reported earlier^’^ from the research collaboration work by NASA and the Nation- 
al Research Council (NRC) of Canada. Further, preliminary analysis of test data from the 
2013 full scale turbofan engine ice crystal icing test^ conducted in the NASA Glenn Propul- 
sion Systems Laboratory (PSL) has also suggested that (1) both types of ice formation oc- 
curred during the test, and (2) the model has captured some important qualitative trend of 
turning on (or off) the ice crystal ice formation process in the tested engine low pressure 
compressor (LPC) targeted area under different icing conditions that ultimately would lead 
to (or suppress) an engine core roll back (RB) event. 
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Nomenclature 

relative heat factor, dimensionless 
modified relative heat factor, dimensionless 

specific ratio of ice particle kinetic heating to latent heat absorbed from melting, dimensionless 
specific heat of air, cal/g K 

specific heat of water at the surface temperature, cal/g K 

convective heat-transfer coefficient, cal/sm^K 

gas-phase mass-transfer coefficient, g/sm^ 

ice water content, g/m^ 

liquid water content of melted ice, g/ m^ 

ice crystal melt ratio, = LWCt/IWCp dimensionless 

local Mach number, dimensionless 

ice crystal median mass diameter, pm 

water droplet median volume diameter, pm 

mass flux of water evaporated per unit time, Ibm/fl^ s 

mass flux of ice/liquid- water particle impinged per unit time, Ibm/fl^ s 


Principal Research Scientist, 21000 Brookpark Rd., MS 11-2, AIAA Associate Fellow. 

^ Research Scientist, Icing Branch, 21000 Brookpark Rd., MS 11-2, AIAA Senior Member. 
^ Aerospace Engineer, Icing Branch, 21000 Brookpark Rd., MS 11-2, AIAA Member. 


I 

American Institute of Aeronautics and Astronautics 


TYL 

melt, imp 

mo 

no 

P 

Pv.WOO 

Pv,ws 

Qconv 

^evap 

^freeze 

Q kinetic 

^melt 

RH 


RH hyg 

RTDl 

SH 

Tf 

Tice 

Tmelt 

T 

^ s 

'^wall 

Twb 

V 


VWQ 
j3o 
<t> 

Af_A„ 


gam 


X 

Y 

Subscripts 


i 

t 

s 

0 

GO 


mass flux of water-melt impinged per unit time, Ibm/fl^ s 

surface melting fraction at stagnation, dimensionless 
surface freezing fraction at stagnation, dimensionless 
pressure, nW 

the saturation vapor pressure of water in atmosphere, nW 

the saturation vapor pressure of water over icing surface, nW 

surface heat loss due to convection, Btu/hr ft^ 

surface heat loss from evaporation, Btu/hr ft^ 

surface heat gain from release of latent heat of fusion, Btu/hr ft^ 

surface heat gain from kinetic energy of ice crystals and water drops, Btu/hr ft^ 

surface heat loss from release of latent heat of fusion, Btu/hr ft^ 

relative humidity, % 

relative humidity measured when cloud on, % 

total air temperature measured in the test section when cloud on, °C 

specific humidity, g of water vapor/kg of dry air 

freezing temperature of bulk water, °C 

ice crystal temperature, °C 

temperature of water melt, °C 

icing surface temperature (=7}), °C 

wall temperature, °C 

air wet bulb temperature, °C 

air static temperature, °C 

air velocity, kt 

vapor water content gained (i.e. = VWCt~ VWCt - [pair SH]^ - [p^/^SH]/), g/ m^ 

collection efficiency at stagnation line, dimensionless 

ice/water mass energy transfer parameter, °C 

latent heat of freezing and evaporation, cal/g 

air/vapor energy transfer parameter, °C 

ice crystal cloud kinetic heating parameter, dimensionless 

ice crystal cloud melting profile, dimensionless 


inlet, initial or dry condition 
target, terminal or wet condition 
at the surface 
stagnation value 
local static condition 


Superscripts 

* = threshold value 

~ = threshold value 


I. Introduction 

N umerous turbofan engine power-loss events have been observed (some as recent as 2013^) on both commuter 
and large transport aircraft near large deep convective weather systems at high altitudes'^. The deep convective 
updraft core can bring large amounts of condensed water from the surface to the cloud top. In theory the pre- 
dicted maximum ice water content (IJVC) above the freezing level in these clouds could reach up to about 9 g/no? at 
around 30,000 ft. In addition, from more recent measurement data available on microphysical properties of deep 
convective clouds it appears that fully glaciated ice crystals have been dominated in the deep convection above the 
freezing level and the ice mass has been concentrated at sub-200 pm sizes, and very little super-cooled liquid water 
content (LWC) has been observed. 

^ Aviation Safety Network News, “Russia investigates serious engine ieing ineident on Boeing 747 -800,” the link to the news 
website is: http://news.aviation-safetv.net/2013/08/20/russia-investigates-serious-engine-ieing-ineident-on-boeing-747-800/20 , 
reported on August 2013 
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Mason et al.^ have proposed a hypothesis for possible ice accretion due to ice crystal ingestion based on the in- 
strumented flight-testing data collected in commuter and transport aircraft engine events. The hypothesis basically 
recognized that ingested ice crystals could begin to melt causing mixed-phase conditions within the compression 
system passage of the engine. This mixed-phase water mass would impinge and accrete on unheated surfaces within 
the engine where the local air temperature was significantly higher than freezing before entering the cloud. This ice 
crystal icing inside the engine compression system has been shown to create flow blockage and hardware damage, 
as shown in Ref 3, such that the engine could no longer function properly. The hypothesis specifically stated that 
liquid water is a necessary condition locally for ice accretion to continue. Furthermore, it was believed that this ice 
crystal accretion phenomenon is not unique to commuter or transport aircraft engine type and therefore it was ex- 
pected that certain common physical mechanisms for ice crystal accretion on all jet engines should exist. 

To better understand the ftmdamental aspects of ice crystal accretion physics and possible influence of high con- 
centrations of ice crystals on engine performance and operability, NASA and NRC of Canada have begun a series of 
ftmdamental icing physics studies^’ ^ to examine the physical mechanisms of ice accretion on surfaces exposed to 
ftilly glaciated ice crystal and mixed-phase conditions similar to those believed to be encountered during engine 
power-loss events. These investigations were mainly experimental utilizing the NRC Research Altitude Test Facili- 
ty (RATFac) and a special cascade wind tunnel, see Figures 1-2 , to create semi-realistic icing local flow conditions 
inside a turbofan engine low pressure compressor (LPC) region to (1) observe how ice crystal accretion process 
takes place, (2) identify the onset conditions and the corresponding threshold limits to turn on or turn off the accre- 
tion and finally (3) explore the inter-relation of relevant physical variables associated with engine ice crystal icing 
processes. 



Figure 1. Ice crystal test system with cascade wind tunnel installed in NRC RATfac 
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Figure 2. The cascade wind tunnel showing experimental configuration including airfoil 


Key findings from those early studies included: (1) the air wet bulb temperature, was identified to be an ef- 
fective parameter to represent the air/water vapor mixture (i.e. humid air) temperature for studying ice-crystal icing 
process. In particular, it was shown that higher altitude (i.e. decreasing pressure) with lower humidity inside an en- 
gine with ingested ice crystals enhances the evaporative cooling of the local air flow and the compression system 
passage wall temperatures and reduces the initial wet bulb temperature to a lower asymptote in a time scale of about 
one half to a few minutes into the spray, (2) two distinct types of ice formation were observed on an unheated icing 
surface. When the (total) air wet bulb temperature T^b,o was above freezing (> +2°C) there seemed to exist a corre- 
sponding IWC threshold (/JVC*) in order to have some ice formed on the surface by the accumulation of un-melted 
ice crystals with relatively weak surface bonding. As a result, a large piece of shed ice was frequently observed. 
Furthermore for any given /JVC > IWC* at the prescribed T^b,o, the overall ice size and the leading edge ice growth 
rate (before any big shedding) became larger as IWC increased. Interestingly it was also observed that for the larg- 
est IWC value tested (i.e. 17 g/m^), there was a corresponding wet bulb temperature threshold (T^b/) that when ex- 
ceeded could eliminate the ice formation from the surface. For example Runs 522 & 524 of Ref 2, the T^b,o was 
raised from about 2.7°C to 9.3°C to suppress the ice formation completely as shown in Figure 3. 

However when Ty^b,o was just below freezing ( — 2°C) the ice was formed by accretion of water melt (similar to 



M,„=0.25, Po=6.5 psia, IWC=17 g/m^ T„b,o=2.7°C. M^=0.25, Po=6.5 psia, IWC=17 g/m^only, T„b,o=9.3°C. 


Figure 3. Large ice formation with frequent shedding was observed in Run 522. For Run 524 no ice was 
observed by raising the Ty^b,o to 9.3°C. 

the super-cooled liquid water icing process) with strong adhesion to the surface. No large ice shedding was noticed 
but sizable mass of ice and water (e.g. small glaze ice feathers and surface runback water) was stripped off from the 
ice structure formed in part due to the known effect of ice erosion and surface water splash-off. This loss effect was 
further amplified as local Mach number increased. For example Runs 512 & 504 of Ref 2, the T^b,o was -2°C, and a 
very glazy ice accretion was observed for Run 512 at Moo = 0.25 but a rime-to-glaze transition type of ice accretion 
was observed for Run 504 at Moo = 0.4, as shown in Figure 4. 


^ Airfoil leading edge 



Figure 4. Two different ice accretions were observed for Twb,o=-2°C at local Moo=0.25 & 0.4 (top view). 
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Run 512 



Figure 5 Evolution of total air temperature and relative humidity measurements in the test section. 

The calculated leading edge ice thickness from LEWICE2D without the mass loss was an order of magnitude 
larger than the experimental images indicated. This type of mass loss phenomenon has been reported earlier by 
Kamel Al-KHalil^ in his 2003 experimental study of mixed-phase icing effect on thermal ice protection systems 
(IPS), and lastly (3) there seemed to be a strong thermal interaction between the ingested ice crystals and the initial 
warm local airflow environment. As the IWC cloud entered the test section, the dry air total temperature started 
dropping and the water vapor content was rising to a new asymptotic value, as shown in Figure 5 for Run 512 of 
Ref 2. This observation was significant and implied that the dry conditions, measured prior to the start of an icing 
test, do not accurately represent the conditions near the test article while the cloud is on. 

The focus of the current study is to present the thermodynamic model developed specifically for engine ice crys- 
tal icing and some analysis results from the model will be validated with aforementioned experimental data trend 
from Ref 1 and 2. Finally, a selected number of engine LPC local flow simulation results from COMDES^’^, Glenn- 
HT, LEWICE3D and surface temperature measurement data from Ref 3 will be used to assess the validity and ap- 
plicability of the model to simulate critical ice formation processes related to a full engine icing event. The validat- 
ed model will become the theoretical ground for new altitude scaling method development^ to support future scaling 
tests as well as preparatory experimental work at RATFac and PSL including ice-crystal cloud characterizations, 
new test rig and new test article development, and new imaging methods for ice shape characterization needed in 
altitude scaling tests. 


II. The Thermodynamic Model for Ice Crystal Icing 

This new thermodynamic model is primarily based on the scaling work originally done by Ruff^ and later ex- 
tended by Anderson^^ and Anderson and Tsao^^ for super-cooled liquid water icing in both Appendix C and SLD 
regimes. The expressions for the parameters used will be presented here without much discussion. Therefore, read- 
ers who are interested in the physical descriptions and detailed derivations of these parameters are referred to Ref 
10 and 1 1 and the references given therein. Some modifications are made in the surface energy balance equation to 
account for possible interactions between the local air flow within the compression system flow passage and the 
ingested ice crystals from the engine inlet and they will be described in the discussion below. 

A. Model Assumptions 

Currently the model only considers the engine ice crystal icing phenomenon for a fully glaciated ice crystal 
cloud which is the most prevalent cloud characteristic in or around a large deep convective weather system at high 
altitude"^. Figure 6 illustrates the notional schematics of ice crystal particle path through an engine fan and LPC re- 
gion. When the surface energy balance analysis is applied to a selected stagnation location of a LP compressor 
component, e.g. a geared turbofan engine 2^^ exit guide vane leading edge (i.e. EGV2 LE), the model assumes or 
considers: 

• local air flow is steady relative to accretion process, 

• accretion is quasi-steady, 

• frilly glaciated ice particle cloud initially with single particle size. 
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• no LWC supplement (i.e. Mly glaciated ice crystal cloud at inlet) 

• when ice particles experience some melting, Tfce = Tmdt = Tf (i.e. 0 < MR < 1 ), 

• when MR = 0 (no melting), Ticej < Tice < Tf, 

• when MR - 1 (hilly melted), Tf < T^eit < To , 

• for partially-melted ice particle cloud, LWCt = IWQ * MR, IWCt = [IWQ * (7 - MR)] - VWCgain. 

• on unheated surface {T^^aii could be fairly hot initially), 

• Ts = 7/ when glaze ice accreted or accumulated on surfaces. 

B. The Surface Energy Balance Equation 


Once ice cloud is on, air starts cooling down, SH LWC^, M VD|. { - lpc condition) 

increasing and then asymptotic to a steady state !WCj, MM D | 



Figure 6. The schematics of ice particle path through an engine fan and LPC region. 

As duly noted in Ref 1 and 2, there are two distinct types of ice growth mechanism in ice crystal icing. It is ad- 
vantageous to present the corresponding equation accordingly from applying the mass and energy conservation laws 
on the icing surface. Here only the final governing equations and key expressions of parameters are presented with- 
out much discussion for brevity. These two types of ice growth mechanism are described as follows: 

(1) Freezing dominated process on surface when cooling > heating (i.e. 0 < ^1) 

^evap ^conv ^kinetic ^freeze ( 1 ) 
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Substituting the expressions for various heat flux terms listed above into Eq. (1) it can be rewritten in the follow- 
ing form 
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The surface freezing fraction is formally defined as the ratio of the amount of water that freezes at a given sur- 
face location to the total amount of water that impinges at that location. The key terms in this formulation include (j) 
and Go, which have dimensions of temperature and relate to the ice/water mass energy transfer and air/vapor energy 
transfer, and bo, the relative heat factor, which was first introduced by Tribus, et al.^^ 
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With some algebraic manipulations, Eq. (2) can be rewritten as 

n„ = UT,-T^,.„)-c, 


(3) 


in which denotes the air total wet bulb temperature and is given as 
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One can show from Eq. (3) that both Ty^b,o and MR can affect the accretion process but they do not uniquely de- 
termine the value of ng. The no depends on other important icing surface quantities like he, fo and nioss^ There are 
many combinations of T^b,o and MR that could have the same freezing fraction ng. It is also noted that for a given 
local thermal profile, i.e. fixed T^b,o and MR, the value of freezing fraction ng will decrease as one increases IWC 
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load. Similar behaviors were observed in super-cooled water icing. This strongly suggests that the most important 
parameter affecting the ice crystal accretion process is the freezing fraction. In addition, by letting rin = 0 , a total 
wet bulb temperature threshold is obtained 


r = r 


■ ^0 * ^0 “ 


- • Cq . 


( 4 ) 


It should be noted that some freezing of water melt would occur on surface if <T whereas some melting of 
ice particle would occur on surface >T . 


(2) Melting dominated process on surface when heating > cooling (i.e. 0 < <1) 
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Carrying out similar substitution procedures as described above, the surface melting fraction, formally defined as 
the ratio of the amount of ice that melts at a given surface location to the total amount of ice that impinges at that 
location, at stagnation can be rewritten from Eq. (5) as 
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( 6 ) 


Similar to the freezing dominated process, Eq. (6) shows that both and MR still affect the ice build-up pro- 
cess but they do not uniquely determine the value of mg. Further, for a given T^b,o and MR, as one increases IWC 
load, mo becomes smaller value 0) and the surface ice build-up will increase (both in size and rate). Likewise, 
for a given IWC cloud, as To rises, mg becomes greater 1) and the surface ice formation will diminish. Asymp- 
totically the surface ice growth is proportional to (7 - mg). 

To further understand whether in this regime the ice crystal icing phenomenon could be possibly turned on (i.e. 

< 1) or turned off {mg =1) by simply changing certain key variables, a special case of mn=\ (i.e. fully melted 
condition) is studied with Eq. (6) and it is found that there exists an ice/water mass flux threshold to survive a com- 
plete melting on any given local stagnation surface 
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There are three distinctly different types of physical parameters involved to determine this threshold: 


1 . he, representing the local surface convective heat transfer characteristics. The greater the local con- 
vective heat transfer is the larger the threshold would become. This suggests that with lower he design 
stagnation regions inside the engine LPC the local surface might not have enough thermal margin to 
melt all possible impinging ice crystals for a given engine operation set point. 

2. (Ty^b,o - Ts), which represents the actual melting potential from air total wet bulb temperature to ice 
melting temperature. For a given engine LPC thermal profile, the lower the air pressure and water va- 
por content at higher flight altitude the smaller the air total wet bulb temperature becomes. Conse- 
quently it lessens the available melting potential inside the LPC which means the threshold would be 
lower. Unless the engine will increase its N 1 setting or the aircraft will lower its flight altitude to raise 
its local air temperature, the ice crystal icing is likely to occur on surfaces with lower threshold. 
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3. and Y, which are non-dimensional parameters denoting the ice crystal cloud kinetic heating and its 
melting profile right before impact. The net effect of this term (Y/;^) parametrically shows that the 
threshold will increase its value if the impinging /JVC cloud could experience larger kinetic heating 
due to higher local inviscid flow velocity or the /JVC cloud could have greater melting before impact. 
As a result, ice crystal icing would be more likely to occur at locations where the local inviscid flow 
velocity becomes slower or in engine operation conditions that do not promote more melting of ice 
particles early inside the LPC right before their impact, (e.g. the high altitude test points like DP 908 
and SR-3 from the PSL engine icing test, see Ref 3, 6 and 7) 


Now considering for a given warm air total wet bulb temperature (> 0 ®Q, from Eq. (7), it is shown that there 
exists a IWQ threshold to allow some ice to remain on the surface and it is written as 


IWC* = VWC + 

gam 




( 8 ) 


It suggests that if IWCi <IWC^, no ice could survive the melting but if IWCt > IWO^, some ice build-up would re- 
main on the surface. As for a given IWQ cloud, Eq. (7) also shows that there exists a T^b,o threshold to melt all the 
ice on local surface. This wet bulb temperature threshold is written as 
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= T 




Y K 


( 9 ) 


Similarly it suggests that if Ty^b,o ^ T* , no ice could survive the melting but if T^b,o < T*, some ice build-up would 
remain on the surface. 


III. Model Analysis and Comparison with Existing Experimental Data 

In this section, some model analysis results are compared with aforementioned experimental data trend observed 
from Ref 1-3. The results are further divided into two parts: 

A.NRC RATfac Test Data - Local Flow Simulation 

The Runs 506, 522, 524 and 530 were conducted on March 9, 10 and 11, 2011 separately with condition listed 
below in Table 1. For Run 506, the initial T^b,o was set to be -2 ®C. During the 6 minutes of IWC only spray some 
ice crystal accretion was formed on the wedge model mid-span leading edge surface with no large ice shedding ob- 
served. The thermocouple reading from the leading edge location was fluctuating around (-0.8, 0.2) ®C. Applying 
the Run 506 test condition to Eq. (4) it was shown that with T^b,o ^ -2 °C, MR = 0.1, and estimated nioss ^ 0.92, the 
total wet bulb temperature threshold is about 

f = -0.8°C. (10) 

Since the initial T^b,o set value is less than the threshold, some freezing of water melt would occur on surface. 


Table 1 Run 506, 522, 524 and 530 test conditions 


Test Conditions 



Run 506 
(03-09-11) 

Run 522 
(03-10-11) 

Run 524 
(03-10-11) 

Run 530 
(03-11-11) 








Ice Water Content 

IWQ 


7.0 

17.0 

17.0 

3.0 

Ice Particle MMD 

MMD 

pm 

115 

115 

115 

115 

Liquid Water Content 

LWQ 


0.0 

0.0 

0.0 

0.0 

Total Temperature 

To 


13 

20 

23.5 

23 

Total Pressure 

Po 

psia 

6.5 

6.5 

6.5 

6.5 

Relative Humidity 

RH 

% 

8.0 

13.0 

27.0 

9.0 

Total Wet Bulb temperature 

Twb,0 


-2.0 

3.0 

9.3 

3.0 

Mach Number 



0.4 

0.25 

0.25 

0.25 

Spray Time 

t 

min 

6 

2 

2 

3 

Angle of Attack 

a 

deg 

-6 

-6 

-6 

0 
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For Run 530 the initial was set to be 3 °C. After 3 minutes of IWC only spray, no ice formation was ob- 
served on the test article surface. The surface thermocouple reading from the wedge model leading edge in time 
average was around 4.0 °C. Applying the Run 530 test condition to Eq. (8) it was shown that with T^b,o = 3 °C, 
y^Cgain = 0.935g/m^, MR = 0.3, and estimated riioss = 0.90, the IWC threshold is approximately 

IWC* ( 11 ) 

Since the IWCt value is less than the threshold, no ice would survive (i.e. mo=l, all melted). However if one let 
the IWCi be greater than the threshold, some ice build-up should remain on the surface. This was exactly the ra- 
tionale behind Run 522. It was a 2 minutes of IWCi only spray with everything the same as Run 530 except IWCf = 
17 g/m^ and a = -6 deg. For Run 522 a large ice formation with frequent ice shedding and growing was observed on 
the wedge airfoil surface at the same T^b,o = 3 ®C, see Fig. 3. 

Further, one could ask if it were possible to turn off this large ice build-up process completely by raising the air 
total wet bulb temperature. This was the design objective for Run 524 in which the T^b,o was raised from 3 to 9.3 
°C. The large ice build-up completely vanished and became surface runback water (see Fig. 3). The thermocouple 
reading for this case in time average was around 7.0 ®C. Applying the Run 524 condition to Eq. (9) it was shown 
that with T^b,o = 9.3 °C, VWCgam = 0.967g/m^, MR = 0.41, and estimated riioss = 0.90, the T^b,o threshold is about 

T*=7.9 C (12) 

Since the Ty^b,o value now is much warmer than the threshold T*, all ice should be melted on the surface accord- 
ing to the model estimates and the test data confirmed that there was no ice remaining on the wedge model surface. 

B. NASA PSL Ice Crystal Icing Engine Test - Full Engine Flow Simulation 

A known AlliedSignal ALF 502-R5 turbofan engine RB condition, FLT 850, was chosen as the baseline for this 
analysis. Due to very limited number of instruments installed in the engine LPC region for the 2013 PSL engine 
icing test, some local flow details of the ALF 502-R5 LFOl engine in FLT 850 condition were extracted from the 
baseline test data DP 443 or were calculated from COMDES/MELT, Glenn-HT and LEWICE3D simulations if the 
needed information was not available. In the analysis the main icing site under consideration is the 2^^ exit guide 
vane stator leading edge (i.e. EGV2 LE). 

The proposed key connection to a full engine RB event is that if ice crystal icing were sustained in EGV2 LE a 
RB event would occur for this engine and the stronger this ice build-up process were in EGV2 LE the quicker the 
RB event would occur. On the contrary, if icing is not allowed in EGV2 LE then no engine RB event would likely 
take place. 

Table 2 displays a number of FLT 850 related test conditions that were selected for the comparison and the test 
conditions are all scaled by the baseline, DP 443 . From the table one clearly can see how the change in IWC, N1 or 
Too could affect the engine RB phenomenon. The PSL test data showed that an engine RB event can be turned off 
completely by reducing the IWC load like DP 574, or by increasing the fan speed N1 like DP 735, or by raising the 
ambient air temperature Too like DP 740. The engine RB event can also be further enhanced by increasing the IWC 
intake like DP 577. It should be noted that the DP 729 is another FLT 850 baseline variant with 57% of the baseline 
IWC value that has experienced a slower called RB event. 

For the first four test conditions in Table 2 Eq. (7) was used to determine the ice mass flux threshold. However 
the actual ice mass flux at EGV2 LE location was approximated from the combined Glenn-HT and LEWICE3D 
simulations. The final results were all scaled by the ice/water mass flux threshold of the baseline, DP 443, for data 
sanitization purpose. 


Table 2 FLT 850 related test conditions in 2013 PSL engine icing test 


Description 

DP 

Alt 



N1 

IWC 

RB (Y, N) 

RB (Y, N) Call Time 

FLT 850 

443 

I 

I 

I 

I 

I 

Y 

1 


574 

I 

I 

I 

I 

0.33 

N 

7.1 


577 

I 

I 

I 

I 

1.23 

Y 

0.9 










FLT 850-1 

729 

I 

I 

I 

I 

0.57 

Y 

1.7 


735 

I 

I 

1 

1.06 

0.57 

N 

6.9 


740 

I 

I 

1.02 

I 

0.57 

N 

3.8 
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Similarly for the last two test conditions in Table 2 Eq. (9) was used to calculate the threshold. To deter- 
mine the local air total wet bulb temperature Ty^b,o the local surface temperature measurement data at EGVl and 
EGV2 were also analyzed to find the estimate. The values were again scaled by the baseline T^b,o threshold. The 
calculated results were tabulated in Table 3. It was noticed that for each RB event test condition there was ice build- 
up at EGV2 LE from the model prediction, and similarly for each non RB event condition no ice was allowed to 
form at EGV2 LE. 


Table 3 The model analysis results 


Data Point 

Location 



Icing (Y, N) 

443 

EGV2 LE 

1.74 

1 

Y 

574 

EGV2 LE 

0.57 

1 

N 

577 

EGV2 LE 

2.14 

1 

Y 

729 

EGV2 LE 

1.00 

1 

Y 






Data Point 

Location 

Twb,0 

T 

Icing (Y, N) 

735 

EGV2 LE 

1.003 

1 

N 

740 

EGV2 LE 

1.006 

1 

N 


As stated in the beginning of this comparison, it was believed and hypothesized by the first author that ice crystal 
icing occurring at EGV2 LE location is directly tied to the formation of an engine RB event. The preliminary results 
from the model analysis seem to agree well with the PSL engine icing test data trend. This would indicate there is a 
strong connection between the local ice formation inside an engine LPC at EGV2 LE and the engine RB phenome- 
non. 


V. Conclusion 

A thermodynamic model is developed to describe possible mechanisms of ice formation on unheated surfaces in- 
side a turbofan engine compression system from frilly glaciated ice crystal clouds. It is shown from the analysis 
that, in general, there could be two distinct types of ice formation: (1) when the surface is dominated by the freezing 
of water melt from frilly or partially melted impinging ice crystals, the ice structure is formed from accretion with 
strong adhesion to the surface, and (2) when the surface is instead dominated by the further melting of ice crystals, 
the ice structure is formed from accumulation of un-melted ice crystals with relatively weak bonding to the surface. 
The model is able to capture important qualitative trends of the fundamental ice-crystal icing phenomenon reported 
early from the research collaboration work by NASA and the National Research Council (NRC) of Canada. 

In addition, preliminary model analyses are performed to compare the model’s estimates with selected data from 
the 2013 full scale ALE 502-R5 LFOl turbofan engine ice crystal icing test conducted in the NASA Glenn PSL. The 
preliminary results showed that the model is able to capture important qualitative trend of turning on (or off) the ice 
crystal ice formation in the engine low pressure compressor EGV2 LE area under different icing conditions that ul- 
timately lead to (or suppress) an engine core roll back (RB) event. Detailed analyses on other interesting test data 
from the PSL test matrix, e.g. the FLT 855 no RB test point and six revenue service RB field event test points desig- 
nated as RS-1 to RS-6 in Ref 3, will be performed to assess the overall validity and effectiveness of the proposed 
model. New altitude scaling methods derived from this model for simulating ice crystal accretion inside a turbofan 
engine will be developed and evaluated in PSL in the coming NASA engine icing test. 
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